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Abstract 

We apply the proposal of Dijkgraaf and Vafa to analyze J\f = 1 gauge theory with 
SO(N ) and Sp(N ) gauge groups with arbitrary tree-level superpotentials using matrix 
model techniques. We derive the planar and leading non-planar contributions to the 
large M SO(M ) and Sp(M) matrix model free energy by applying the technology of 
higher-genus loop equations and by straightforward diagrammatics. The loop equations 
suggest that the MP 2 free energy is given as a derivative of the sphere contribution, 
a relation which we verify diagrammatically. With a refinement of the proposal of 
Dijkgraaf and Vafa for the effective superpotential, we find agreement with held theory 
expectations. 
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1 Introduction 


Recently Dijkgraaf and Vafa have proposed that the exact low-energy super¬ 

potential of certain Af = 1 supersymmetric gauge theories is captured by the large M 
behavior of certain associated M x M matrix models. This is quite remarkable, as 
it reduces the problem of what is, in general, strongly-coupled physics of the confin¬ 
ing phase of pure gauge theory to the zero-dimensional dynamics of a matrix integral. 
Furthermore, the gauge theory quantities are computed from just the planar graphs of 
the matrix theory, nevertheless capturing finite N results in SU(N) gauge theory. 

The conjecture was initially tested for Af = 2 SU(N) gauge theories softly broken 
to Af = 1 by a tree-level superpotential for the adjoint chiral superfields [llj and for the 
J\f = 1* deformation of Af = 4 SU(N) SYM |J, [|, [5|. The conjecture has since been 


extended to a number of other cases 


22|,^3[ , |24l,|25l , ^6H and has been derived from SU(N) gauge theory p7j|28|,^9|. 

In this work we use matrix techniques to analyze Af = 1 gauge theory with SO(N) 
and Sp(N) gauge groups. By a careful consideration of the planar and leading non- 
planar corrections to the large M SO(M ) and Sp(M) matrix models, we derive the 
matrix model free energy. We do this both by applying the technology of higher-genus 

and by straightforward diagrammatics (see e.g. 102,[33||). 


loop equations of 00 


The outline of the paper is as follows. In section we discuss general features 
of the four-dimensional gauge theories with Af = 2 supersymmetry softly broken to 
Af = 1. We also suggest a result for the superpotential of the J\f = 1* theory with 
gauge group SO(2N) which is based on a generalization of a derivation of Dorey 


for SU(N) gauge group. In section |3j, we discuss the geometric engineering of the softly 
broken Af = 2 gauge theories by wrapping D5-branes and 05-planes on compact cycles 
of generalized conifolds. 

The corresponding matrix models are introduced and solved in the large M limit 
in section |]. As for SU(M), we find that the loop equation for the resolvent of the 
matrix model describes a Riemann surface which is identified with a factorization of 
the spectral curve of the Af = 2 gauge theory. The large M solution of the matrix 
models computes the special geometry of Type IIB string theory on the associated 
Calabi-Yau manifold. 

In section [|, we discuss the application of the higher-genus loop equations to the 
computation of the RP 2 contribution to the free energy. The loop equations take the 
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form of integral equations which give recursion relations between the contributions to 
the resolvent at each genus. They suggest a very simple solution for the MP 2 contribu¬ 
tion in terms of the sphere contribution. In fact, the one crosscap contribution to the 
resolvent uq satisfies 

_i du 0 . . 

Ul ~ ±q dSp' ^ ^ 

where Up is the contribution to the resolvent from the sphere. We verify this relationship 
by explicitly enumerating several types of diagrams. We find that the contribution to 
the free energy T\ from MP 2 and from S 2 are related by 

t _u dT{) (1 o) 

Tl = (L2) 

where Sp is half of the’t Hooft coupling for the SO/Sp component of the matrix group. 
We determine the proportionality constant q from the diagrammatics to be q = 

Our results suggest a refinement of the proposal of Dijkgraaf and Vafa for the 
effective superpotential in the case of SO and Sp gauge groups. We find that 

B7- 

Weff = Qd 5 + Q 05 Go ~ 2m r S, (1.3) 


where Qo 5 is the total charge of D5-branes, Q 05 is the total charge of 05-planes, JP 0 
is the contribution to the matrix model free energy from diagrams with the topology 
of a sphere and Go is proportional to T \, the contribution to the free energy from MP 2 
diagrams. We use (|1.3|) to obtain results consistent with gauge theory expectations. 
In particular, the matrix model is consistent with the requirement that there is a 
degeneracy of the massive vacua of the gauge theory given by h, the dual Coxeter 
number of the gauge group. I 11 the case of the AT = 1* SO( 2 N) theory, which we 
discuss in section |(n^, we find that the critical value of the superpotential exactly 
matches the result obtained from our gauge theory arguments in section |2j. We end 
with a discussion of our results and point out several areas for future development. 
Various supporting technical calculations are contained in appendices. 

I 11 the course of this work, two papers on matrix models with SO/Sp groups have 
appeared. In [HJ , aspects of the geometric engineering of the gauge theories, as well the 
leading order in M computation of the free energy of quartic orthogonal and symplectic 
matrix ensembles are discussed. We use a different basis of matrices and we account 
for the appearance of diagrams involving pairs of twisted propagators that are not 
including in the oriented theory. More recently, while this manuscript was in a final 


4 







stage, |J5| appeared. These authors discuss a perturbative derivation of the matrix 
model along the lines of [^j, including a discussion of MP 2 corrections. Their results 
confirm aspects of the refinement (|1.3|) that we found was necessary for the computation 
of the gauge theory effective superpotential. 


2 Results from N = 1 Gauge Theories 


In this section we review some results about Af = 2 and Af = 1 supersymmetric gauge 
theories, specifically the spectral curves and how they factorize when Af = 2 is softly 
broken to Af = 1. We focus on the case of SO and Sp gauge groups (SU was treated 

in 


36| , which this discussion follows). 


2.1 Af = 2 Softly Broken to Af = 1 


As is well known [|37|, |38|[]. the moduli space of Af = 2 gauge theories is governed by 
a “spectral curve”, the periods of which give the masses of BPS objects in the theory 

these spectral curves were found for 


(W-bosons, monopoles and dyons). In ||4T|,|42|, 
the SO/Sp gauge groups. For a rank-r gauge theory, the spectral curve is a genus r 
hyperclliptic curve, 

y 2 = PW 2 O, {</>;}), (2.1) 


where P 2 r +2 is a polynomial of degree 2r + 2 in the x that also depends on the moduli 
0*. The SO and Sp spectral curves can also be written as a genus 2r — 1 curve, 


y 2 = Plri-J- 2 - )>(}), 


( 2 . 2 ) 


which is therefore symmetric under the Z 2 action x 1 —► — x and is a double cover of the 
genus N curve O via this map. 

I 11 Af = 1 language, the Af = 2 vector multiplet of the Yang-Mills theory is decom¬ 
posed into an adjoint chiral superfield $ and an Af = 1 vector superfield V. Af = 2 
supersymmetry can be broken to Af = 1 by an appropriate gauge-invariant superpo¬ 
tential term for $. Because the trace of odd powers of matrices in the Lie algebra of 
SO(N)/Sp(N) vanishes, in contrast to the U(N) case discussed in the super po¬ 

tential deformation for SO(N)/Sp(N) only includes polynomial terms of even degree: 

1 See for reviews. 
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(2.3) 


wwm = J2 fN* 2 *)- 

fc=i z 

A superpotential kbtree of order 2 n + 2 breaks the gauge symmetry down to a direct 
product of n + 1 subgroups, e.g.\ 

SO(N) -> SO(N 0 ) x U(N{) x ... x U(N n ), (2.4) 


where iV = iV 0 + 2A7 H-b 2iV n . 

The 17(1) factors in this theory decouple in the IR. In the supersymmetric vacua 
of the J\f = 1 theory r — n mutually local monopoles simultaneously become massless 
and condense, leading to confinement of the gauge theory |36|. The condition that 
r — n mutually local monopoles become massless leads to a “factorization locus” in the 
moduli space of the spectral curve, where r — n of the (non-intersecting) cycles of the 
spectral curve are simultaneously pinching off to zero-volume. 


Imposing this condition is therefore equivalent to the factorization |44 


2 n 


y 2 = n ^ 2 ~ p * 2 ) 2 ri (- 2 ~ 

j= i 


(2.5) 


2=1 


where pi ^ p 3 , q t ^ qj for i ^ j. On this locus we then obtain the “reduced spectral 


curve 


V 


2 n 


( 2 . 6 ) 


y 2 = IR 2 _ ?2 )> 

3 =1 

which is a genus 2n — 1 curve. This curve parameterizes the J\f = 2 vacua that are not 


lifted by the deformation to J\f = 1 ( |2.3| ) . Notice that the factorized curves now have a 
similar form for SO and Sp, and the curve is still invariant under x ^ —x (this implies 
that the branch points come in pairs: (—(/,;, Qij)- This reduced spectral curve will be 
derived from string theory in the following section by taking an orientifold action on 
the configuration of D-branes on a generalized conifold that engineers this J\f = 1 gauge 
theory. 

The low-energy effective super potential of these gauge theory can be obtained from 
the reduced spectral curve as discussed by It will take the form 


W es = 


E 


NilL — 27T iTiSi 


(2.7) 
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where 2mSi are the periods of the meromorphic 1-form ydx around the A-cycles of 
the spectral curve, fl, the corresponding periods around the B-cycles, and iVj is 


f N t SU(Ni), 
i f-1 SO(Ni), 
l^ + l Sp(Ni). 


( 2 . 8 ) 


By contrast, we will find that the shift IVj i—► IVj emerges in the matrix model by 
considering the first subleading corrections to the large M expansion, coming from 
Feynman diagrams of topology MP 2 . 

Recently the gaugino effective superpotential has been perturbatively derived from 
SU(N) gauge theory |27|, |2B|, . It is interesting to note that similar arguments^] to 

those of |28| can be used to argue that only diagrams with at most one boundary (if 
quark flavors are present) or crosscap will contribute to the gauge theory superpoten¬ 
tial @M@. 


2.2 The Af = 1* Theories 


The Af = 1* theories arise as deformations of Af = 4 Yang-Mills theory by mass terms 
for the three adjoint Af = 1 chiral fields. The total superpotential is 

w = ti U 1 [$ 2 ,$ 3 ] + ) ( 2 - 9 ) 


2=1 


and the F-flatness conditions can be written as 




( 2 . 10 ) 


Supersymmetric vacua are then obtained by embedding £77(2) representations into the 
gauge group G. In particular, for G = SU(N), the embeddings are classified by the 
divisors d of N, leading to d massive vacua |46| . 

Dorey |34|] (see also |47|]), following the approach of (|8|, compactihed the theory 


with SU(N) gauge group on a circle of radius R. The degrees of freedom of the effective 
2 + 1-dimensional theory are r = rank(G) complex Abelian scalar fields X a which are 
composed of the Wilson lines and the scalars dual to the massless photons of the theory. 
The moduli space of the theory is 


M = E r T /W G , 


( 2 , 11 ) 


2 We thank Jaume Gomis and Jongwon Park for discussions on this issue. The same observation 

about crosscap contributions was made in pfj and more recently for boundaries in |ph|. 
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where E r is the elliptic curve parameterized by each X a and Wq is the Weyl group of 
G. 


By several arguments, including the relationship between the elliptic Calogero- 
Moser systems and the J\f = 2 theories which can be softly broken to the Af = 1* 
theory, Dorey found that the superpotential of the 2 + 1-dimensional theory took the 
form 

W = cJ2p(X a ~X b ), (2.12) 

a>b 

where p(z) is the Weierstrass function. Dorey argued that the coefficient c was in¬ 
dependent of the radius R, so that critical values of (|2.12|) (which depend on the 
modular parameter r and not the X a ) could be evaluated in the vacua of the theory 
and extrapolated directly to the R —> oo limit. 

It is interesting to ask what the generalization^ of (|2.12 ) is to arbitrary gauge 
groups G. The modular properties of the superpotential that were crucial to Dorey’s 
argument must still be preserved, so the superpotential should remain a sum of Weier¬ 
strass functions. An obvious guess for the argument of these functions is to replace 
X a — X b by the sum over the positive roots Y^ a >o a • A". The integrable systems ap¬ 
proach |5(],[5]],|46|,[52|] to Af = 2 theories is a promising route to this result. In fact, 
D’Hoker and Phong |53],[54],|55|] have determined the integrable systems that govern a 


large class of Af = 2 theories with gauge group G. An application of the techniques 
of [[K||3]§ to the soft breaking of these theories to Af = 1* suggests that the correct 
superpotential for gauge group G is 


W = c 


{Qi>0} 


p(o L ■ X) 


E 

{ag>0} 


Pu(as)(,^s ■ A) 


(2.13) 


where atL,s are the long and short positive roots of the Lie algebra of G, respectively, 
and p v {z) are the twisted Weierstrass functions 

v — 1 

(2.14) 


pv{%) y ' p t 2 u a 


( 7=0 


defined in |53[. For non-simply laced groups, roots of only two different lengths appear: 
u(a) = 1 for all long roots, u(a) = 2 for all short roots of b n , c n f 4 , while v(a) = 3 for 
the short roots of g 2 - 

3 After an initial version of this paper appeared, we became aware of earlier work by Kumar and 
Troost Q, where they also suggest the formula ( 2.13 ) and give many more arguments for its validity. 
We thank P. Kumar for making us aware of this. 















For SO(2N), since it is simply-laced, twisted Weierstrass functions do not appear, 
and we obtain 

W = cJ2H X a ~ X b) + p( X a + X b)\ ■ (2.15) 

a>b 

Following 0, we can evaluate this in the k th confining vacuum to find (up to an 
additive constant) 

W ~ E 2 ((r + k)/(2N~2)), (2.16) 

where E 2 (t) is the second regularized Eisenstein series. It is tempting to conjecture 
that the result for arbitrary G will take this form with the obvious substitution of h, 
the dual Coxeter number of G for 2 N — 2, but this remains to be verified. 


3 Calabi-Yau Geometry 


We now review the string theoretic engineering of a softly broken Af = 2 gauge theory 
with SO/Sp gauge group |4^,^6|,|56|,[1]. We consider type IIB string theory compactihed 
on the non-compact A\ fibration 


u 2 + v 2 + w 2 + W'(x ) 2 = 0, 


(3.1) 


where W(x) is a degree n + 1 polynomial, which will later be related to the tree level 
superpotential. This fibration has singularities at the critical points of W{x). In the 
neighborhood of those singularities, we can introduce the coordinate x' = W'{x). Then 
it is easy to see that the singularities are all conifold singularities. 

This generalized conifold can be de-singularized in two ways: it can be resolved or 
it can be deformed. The resolution is given by the surface 

u + iv w + iW'(x) \ f \i \ . 

= 0 (3.2) 

—w + iW'(x) u — iv I \ A 2 J 

in C 4 x P 1 . In this geometry each singular point is replaced by a P 1 . These P 4, s 
are disjoint, holomorphic, have the same volume and are homologically equivalent. 
The latter property can be seen by making use of the fibration structure away from 
W'{x) = 0. This Ai fibration over the x plane induces a fibration of non-holomorphic 
S 2, s over the x plane. This S 2 cannot shrink to zero size as one approaches a critical 
point of W in the x plane, but it becomes the holomorphic P 1 of the resolution. 

We can now construct a softly broken J\f = 2 U(N) gauge theory with tree level 
superpotential W (x) by wrapping N D5-branes around the S 2 . This is an UV definition 
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of the theory. A classical supersymmetric vacuum is obtained by minimizing the volume 
of the D5-branes. This amounts to distributing a collection of N t D5-branes over the 
n minimal-volume holomorphic P^s at the critical points of W. The U(N ) gauge 
symmetry is then spontaneously broken to U{N\) x • • • x U(N n _i). 

We now want to consider an orientifold of this theory[]. Since we started with a type 
11B theory on a Calabi Yau, we have to combine the worldsheet orientation reversal 
with a holomorphic involution of the Calabi-Yau (an anti-holomorphic involution would 
be appropriate for the IIA theory). Furthermore we want to fix one of the P 1 ’s and 
act freely on the rest of the Calabi Yau geometry. This can be done if W(x ) is an 
even polynomial of order 2 n. In terms of the fibration structure of the Calabi-Yau this 
means that the critical points of W'^x) come in pairs (— Xi,Xi) and one critical point 
is fixed at x 0 = 0. Then 


(u, v, w, x , Ai, A 2 ) ( -u , -v, -w, -x, Ai, A 2 ) (3.3) 

is a holomorphic involution of the geometry (|3.2|) , which leaves only the P 1 at u — v — 
w = x = 0 fixed. In the string theory this means that there is an 05-plane wrapping 
this P 1 in the Calabi-Yau geometry. 

There are essentially two choices of 05-plane with which we can wrap the fixed P 1 . 
They are distinguished by a different choice of worldsheet action and carry RR 5-form 
charge of ±1 (the RR charge of an Op^-plane is ±2 p-5 in conventions where we count 
the charge of N/2 D-branes but not their N/2 images). The orientifold contribution 
to the RR charge of objects wrapping the P 1 will cause a shift in the coefficient N 0 in 
the flux-generated superpotential on the deformed Calabi-Yau geometry, as explained 
below. 

Now we can construct a softly broken J\f = 2 SO(N)/Sp(N/ 2) gauge theory with 
tree level superpotential W/x) by wrapping N D5-branes around the S 2 and then 
performing the orientifold. The gauge symmetry is again broken SO(N ) i—> SO(N 0 ) x 
U(N 1 ) x ■ • • x U(N n _ ± ) or Sp(N/ 2) ^ Sp(N 0 / 2) x U(N X ) x • ■ • x [/(A^) respectively 
with N = N 0 + 2 Ni + • • • + 2A n _ 1 . 

If we flow this ultraviolet theory to the infrared, there will be a confinement transi¬ 
tion. In string theory this is described by a geometric transition in which the resolved 
conifold geometry with wrapped D5-branes and 05-planes is replaced by a deformed 

l,|5S|l, while the discussion of is more 

closely related to the B-model which is our interest here. 


4 Orientifolds were discussed in the A-model in 57, 56 
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conifold geometry 


IFTTl 


u 2 + v 2 + w 2 + W'(x ) 2 — /(x) = 0, 


(3.4) 


where f[x) is an even polynomial of degree 2 n — 2. Such a polynomial represents 
the most general normalizable deformation of the singular conifold that still respects 
the holomorphic involution ( |3.3| ). For a reasonably small f(x), each critical point of 
W'(x) is replaced by two simple zeros of W'(x) 2 — fix). This means that each P] is 
replaced by a 3-sphere A t with 3-form RR-flux H through it, equal to the amount of 
D5-brane and 05-plane charge on the Pi. The orientifold acts on one 3-sphere A 0 as 
the antipodal map, while the other 3-spheres are mapped to each other in pairs A j and 
A-i. Note that there is no orientifold fixed plane anymore. 

The coefficients in f{x) are normalizable modes that are localized close to the tip 
of the conifold. The coefficients in f{x) are determined by the periods 


Si = - / a 

2 m Ja< 


(3.5) 


These periods Si can be interpreted as the gaugino condensates of the gauge theory. 
There are non-compact 3-cycles Bi that are dual to the Ai. The periods of the B-cycles 


are 


aPo 

dSi 



(3.6) 


where JF 0 is the prepotential. One needs to introduce a cutoff in order to make these 
periods finite. 

The flux through the cycles Ai is determined in terms of the RR-charges of the 
D-brane and O-plane configuration 


N 0 ±2 = 

f H ’ 


Ja 0 

N t = 

1 H ’ 


JAi 


(3.7) 


and the flux through the cycles Bi is given in terms of the coupling constants 


H. 

J Bi 


(3.8) 


Since there is no orientifold fixed plane, there are no contributions to the effective 
superpotential for the gaugino condensate from unoriented closed strings [Q. It is 
then given by the flux superpotential p2|, 


WeffiSi) = H Aft, 


(3.9) 
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where the integral is taken only over half of the covering space of the orientifold. Using 
the expressions for the periods and the fluxes and taking into account the orientifold 
projection, we get 


W'eff(S) 



£Vo 

as, 


«>0 


dfo 

dS r 


l ToSo 




i>0 


(3.10) 


This result could also have been computed on the open string side before the tran¬ 
sition. On the open string side there is no flux through any 3-cycles, so there is no 
contribution to the superpotential due to closed oriented strings. But there are two 
kinds of other contributions to the effective superpotential: the open string contribu¬ 
tions (disk diagrams) and the contributions due to closed unoriented strings at the 
orientifold fixed plane (MP 2 diagrams). The contribution due to the open strings is the 
equal to one half that of the theory without the orientifold, he., it is 

(am 

u i> 0 1 i> 0 


The contribution due to the unoriented closed strings then must be 

WeK(Si) = W ef f(Si) - WZ ff (Si) = ±^. (3.12) 

We will confirm this result in our matrix model computation. 


4 The Classical Loop Equation 


We first consider the saddle point evaluation of the one matrix integral for SO(M) or 
Sp(M) matrices. Our discussion is analogous to that of ]6(||T| and consists of obtaining 
a loop equation for the resolvent. In the next section, we will formulate a systematic 
method to obtain the g s corrections to the classical solution. 

The partition function for the model with one matrix $ in the adjoint representation 
of the Lie algebra of G = SO(M ) or Sp(M ) is 


z= vkdl d * exp {-^ wm )- (41) 

In Appendix [A.l| , we collect results that are useful for SO/Sp groups, but here we shall 
discuss only the SO(2M) group in detail. 
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In the eigenvalue basis, the integral over an SO(2M) matrix is given by 


,. M 


(4.2) 


i=l i<j 

The effective action for the gas of eigenvalues is given by 


5(A) = - e ln < A . 2 - y 2 ) 2 +- E 

• ,. 9 s 

i<j i 


(4.3) 


Note that W is now a polynomial of order 2n with only even powers. This is because 
the trace of an antisymmetric matrix vanishes. In principle 1T( < 1 ) ) could also contain 
the Pfaffian, but we will omit this case. 

This action gives rise to the classical equations of motion 


E - ~W\K) =0. 


A? - Af g. 


(4.4) 


It is useful to define the resolvent 


w o(a0 = .9*Tr = g s ^ 


x — $ 


(4.5) 


then, by multiplying the equations of motion by 2 2 ^\ 2 and summing over i, we obtain 

X A i 

an equation for Uq(x) exactly as for SU(M ): 


u 0 (x) -g, 


2 / (a) 


x 


Uq{x) ) + f(x) - 2lo 0 (x)W’(x) = 0, 


where 


/ 0 ) =9s^2 


2\ i W\\ i ) - 2xW'(x) 


A 2 — x 2 

l L 

is a polynomial of order 2 n — 2 with only even powers, he., it has n coefficients. 


(4.6) 


(4.7) 


In the small g s limit, Q4.6|) reduces to 

u 0 (x) 2 + f(x) - 2 uj 0 (x)W > ( x) = 0, 


or 


where 


V 1 - W'(xY + f(x) = 0, 


y(x) = uo{x) - W‘{x). 


(4.8) 

(4.9) 
(4.10) 
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The force equation is then 


9 m dS 
2y(A) = -9sq£, 


(4.11) 


where the factor of 2 comes from the fact that the force is acting on an eigenvalue and 
its image. This is the same equation as for SU(2M ), with the only difference being 
that the polynomials W and / have only even powers. This matches the expected 
result from the orientifold procedure in string theory. 

The equation for the resolvent can be solved using ( [4.9|) , yielding a formal solu¬ 
tion fib] 

cu 0 (x) = W\x) - V W'[x ) 2 - f(x). (4.12) 


The resolvent is thus expressed in terms of the n unknown coefficients that appear in 
the polynomial f(x ) defined in ( [I.7| ). From the form of the solution, it is clear that the 
resolvent has branch cuts among which the eigenvalues of the matrix are distributed. 
In the large M limit, we thus get a distribution of eigenvalues, with the eigenvalue 
density given by p{ A) 


u; o 0r) = 2 

Jo 

which implies that 
1 


xp(X)dX 


x* 


A 2 


p(X)dX 


x 


— H-— 

A x T A 


p(X)dX 
x — A 


p(A) = — (oj 0 (A + iO) - u> a (\- sO)) = ^—(yj + *0) - y(\ - iO)). 


(4.13) 


(4,14) 


The filling fractions are then given by 

-L f 

Ani 
1 


^n = 


-4o 


Si = 


2 ni 


y(x)dx, 
y(x)dx , i > 0 


(4.15) 


Note that we only integrate around half of the cycle Aq because of the orientifold 


projection. At the classical level, one can see from ( 4.11 ) that y(x) is the force acting 
on an eigenvalue. Now, the variation of the free energy JF 0 of the matrix model caused 
by a changing the number of eigenvalues on the I th cut is then the line integral of the 
force over the non compact Bi cycle of the Riemann surface (4.9) 

dF 0 


dS, 


= / y(x)dx. 


(4.16) 


This is the differential equation that determines JF 0 , be., the leading contribution to 
the free energy. 
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For SO(2M+l) and Sp(M), one can easily see that JF 0 and the Riemann surface are 
the same as in the case of SO(2M). In the next section we will determine the leading 
contribution from unoriented diagrams to the free energy, which is a subleading term 
in the g s expansion of the free energy. 

5 g s Corrections and Loop Equations 

The partition function of the SO/Sp matrix model is 

Z = e 7, r = J dte-A® '1*1 (5.1) 

where the overall coupling constant g s can be thought of as the string coupling and the 
action is 

OO 

W(%) = J2 (5.2) 

3 =1 J 

Dijkgraaf and Vafa IBS conjectured that the exact superpotential of the gauge 
theory with the tree level superpotential W (<h) is given by the perturbative expansion 
of the matrix integral ( |5.1|) around one classical vacuum (saddle point). Such a classical 
vacuum is given by a distribution of the eigenvalues of over the critical points {x*} 
of the superpotential W(x). We denote the number of eigenvalues at the critical point 
Xi by Mi and define the parameters 

So = 9s-^~, Si = g s Mi. (5.3) 

The perturbative expansion around such a classical vacuum can be visualized in terms 
of fat graphs, where edges of a ribbon correspond to Chan-Paton factors. For each 
Chan-Paton factor we have to choose a critical point x t , on which it sits and a loop 
of such a Chan-Paton factor gives a contribution of M* = Si/g s . From the overall 
normalization of the action, it is clear that each vertex of the diagram contributes a 
factor of 1 / g s and each propagator contributes g s . Thus the overall power of g s counts 
the Euler character of the fat graph. The superpotential is then given only by the 
contributions from planar and MP 2 diagrams. These diagrams have a very simple g s 
dependence, but the Si dependence can actually be quite complicated, as we will see. 

If one did the full matrix integral, there would be a sum over all saddle points 
and the Si dependence would be lost. However, since we are interested in only in a 
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perturbative expansion around a classical vacuum, the Si dependence is nontrivial and 
will describe how the effective superpotential depends on the gaugino condensates. 

In a recent paper [67j, it has been shown that the in a vacuum where the gauge 
symmetry is broken to a subgroup (say, a product of U(Ni ) factors), the off diagonal 
components of the matrix <f> do not correspond to propagating degrees of freedom, 
and that these should be properly interpreted as the Faddeev-Popov ghosts that are 
necessarily included because of the gauge fixing involved in doing the matrix model. 
Thus, for computing Feynman diagrams in the matrix model, we have to include terms 
in the Lagrangian that belong to the ghost sector as well. But the loop equations, as 
we shall see, correspond to Ward identities in the matrix model. They arise because 
of the invariance of the matrix integral under an arbitrary reparametrization of $ that 
respects the SO/Sp symmetry of the Lagrangian. If we take into account the variation 
of the measure as well, then this symmetry leads to the loop equations. Thus, we do 
not expect the ghosts to be relevant for the discussion in this section. 

In the SO/Sp case we expand the matrix model partition function in a systematic 
expansion in g s . The coefficients of the terms in the expansion are the contributions 
coming from the Feynman graphs that can be drawn on a surface of Euler character 
X = 2 — 2 g — c where g denotes the genus, and c denotes the number of cross-caps. 
We mentioned earlier that each loop in a Feynman diagram contributes a factor M. 
In order to see this, consider the propagator for the SO(M) matrix model. It has a 
group theoretic factor 

(®ij$ki) ~ -(fiikdji ~ Siibjk)- (5.4) 

Thus, each loop in a Feynman diagram contributes a factor of M t (the number of 
Chan-Paton factors on the ?’th critical point). 


5.1 The Resolvent 


We shall now introduce the general technique of loop equations, which is an iterative 
procedure to calculate the higher order (in g s ) corrections to the partition function. 
Central to this procedure is the loop operator defined as 


d _ 2 j d 

dV [X) ~ ^ aW +1 dg,' 
j =i J 


(5.5) 
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The resolvent, which is the generating functional for the single trace correlation func¬ 
tions of the matrix model is defined as 


( , 6) 

' 1 k =0 

Using the identity 

-{2k)^-T = g s {Tx $ 2fc ), (5.7) 

dgk 

we can express the resolvent as 

u(x) = + -, (5.8) 

dV x 

where we used S = = g s M. We are using the variables g s and S, since we are 

working in the small g s limit with S fixed. As mentioned before, the free energy has 
an expansion in g s of the form 


J r =E». 2 " +0 ^ 

9,c 


(5.9) 


We will be interested in calculating the first two terms in this expansion, which are 
the contributions from diagrams with the topology of S 2 and MP 2 . The resolvent has 
a similar expansion 

u(x) = ^g 2 s 9+c uj gfi {x). (5.10) 

9,c 

The asymptotic behavior at infinity of the cn 9)C is clear from the definition of uj(x) 


^o,o(^) =- + 0 {x 2 ), 
x 

UgA x ) = ^(x 2 ), 2g + c > 0. 


(5.11) 


Using this fact and the existence of the genus expansion, we can write 

d S 

U 0 ,o( X ) =-Jy( X )d r 0,0 + 

Ug,c( x ) =-^( x )J 7 g,c 2g + C> 0. 

These equations determine the dependence of T g ^ c on the coupling constants. There 
is still an additive constant which is undetermined, but this is unphysical. In the next 
section we will derive the loop equation, which will provide us with recursion relations 
to calculate cn 9jC as functions of the coupling constants gj. (For the rest of the discussion, 
we denote cno,o by uq and ojq.i by u q.) 
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5.2 The Loop Equation 


In this section we will derive an important recursion relation between the different 
perturbative contributions a; 9iC to the resolvent^. The loop equation can be derived by 
doing a reparametrization of the matrices $ in the matrix integral and observing that 
the integral is invariant under this reparametrization. Let us reparametrize $ by 






oo 2 

‘“'-'EE 


odd 

<y l dsl>'¥ 


~ &2k+l 

* - £ E 


k =0 


x 


2 k+2 


k =0 «=0 


X 


2fc+2 


(5.13) 

(5.14) 


where we only take the odd/even powers of <L / in order to preserve the SO/Sp Lie 
algebra. The Jacobian for this reparametrization is then 


J($') = 1-1 (Tr — 

V ; 2 V x - <L' 


+ — Tr 1 


2 x x — 


The action transforms as 


Tr W($) = Tr W — 


x — <&' 


odd 


= Tr W($') - eTr 


JT / ($ / ) 
x — 


Inserting this into the matrix integral, we get 


d& 


Tr 


x — <&' 


- -Tr 1 


x x — <E>' 


- —Tr W($') 
e 9s v ’ 




9s 


We can now make use of the identity 


x — 


(5.15) 


(5.16) 


(5.17) 




Tr 


x — <3> 


(5.18) 


to get the loop equation 


/ W'($)\ 1 g s g 2 s d 

9‘ \ Tr / = 2 w(l) - 2i“ {x) + T5v 1 - xMx) - (5 ' 19) 


5 In an earlier version of this paper, the last term in ( [5. 15| ) was missed. This was corrected in ps| |. 
While our result for the RP 2 contribution (5.34) is unchanged, we have corrected our derivation here. 
The same result was obtained via different methods in 11681 . 
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We can rewrite the loop equation using 


9s ( Tr 


W'($) 

x — <3> 


9 * < £ mr X,) 


x - A i 


dx' W\x' 
c 2m x — x 


Jtv(x'), 


(5.20) 


where C is a contour that encloses all the eigenvalues of $ but not x. In the small g s 
(large M ) limit of the matrix model, we get a continuous eigenvalue distribution for $ 
and all the eigenvalues are distributed over cuts on the real axis of the x-plane. The 
loop equation now reads 


dx’ W\x’) , 1 / 9 a ( n , 9, - / N / N 

-U)(x ) = -a }(x) — —U)(x) + -^-——(x)uj{ x ). 


d 


2 m x — x 1 


2 x 


2 dV 


(5.21) 


We can now insert the g s expansions for the resolvent and iteratively solve for the u 9iC . 
The zeroth and first order equations are 


dx 1 W'(x') , 1 , , 2 

- ^L w ( X } w 1 (x') = oj 0 (x)oj 1 (x) - ^-CJo(x). 

c 2m x — x 2x 


c 


(5.22) 

(5.23) 


The resolvent that solves the loop equations has to satisfy (|5.11|) which imposes 
constraints on the end points of the cuts in the x-plane. Note that this derivation of 
the loop equation is valid in the saddle point approximation that we are using. 

Equation (|5.23 ) is a linear inhomogenous integral equation for uq. The homoge¬ 
neous equation is solved by a derivative of u; 0 with respect to any parameter which 
specifies the vacuum, i.e., is independent of the coupling constants g 3 . In our case 
there are only the parameters 5), which specify the classical vacuum around which the 


matrix integral is expanded. We will elaborate on this observation in section |5.3.2| for 
the case of a softly broken J\f = 2 theory and we will use this result for the J\f = 1* 
theory. 


5.3 Solution to the Loop Equations 


Let us now solve the the loop equations (|5.22|) first for ojq and then for uq in the case 
of a polynomial potential 

n 

W($) = ^^ : $ 2j . (5.24) 

3=1 2j 


In this section, we closely follow the discussion in [30,31 . 
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5.3.1 Planar Contributions 

In equation ( |5.22D , we deform the integration contour C to encircle infinity, and rewrite 
it as 

i r ap \ai'(p\i 

(5.25) 


-uj 0 {x) 2 = W'(x)u Q (x ) + 


c a 


dx' W' {x')ujq{x') 
2m x — x' 


Assuming that Uq{x) has k cuts in the complex x-plane, we make the ansatz 


loq(x) = W'{x)-M{x) t 


2k 


\ i =i 


(5.26) 


where M(x) is an undetermined analytic function at the moment. Here the end points 
of the cuts, denoted by the x t , are unknown and have to be determined. From the 


discussion in |3l[], it is clear that if we have the maximum number k = 2n — 1 of cuts 
allowed, the function M{x) is a constant. The loop equation determines M in this case 
to be the coupling constant g n . Also, in the SO/Sp case the eigenvalues come in pairs 
and the total number of “independent” cuts is n. There is one cut [—xo,xo] centered 
around zero, and the other cuts come in pairs [x 2 i-i, X 2 i\ and [— X 2 i, —a? 2 i— i]- We shall 
follow these notations in what follows. 

We now demand that the resolvent ujo(x) have the S/x fall off at infinity and thus 
get n constraints 


dk,n _ 


dx' 


X 


,2fc " 1 W’(x ' 


:, fc=l,2, 


2/c2 ”v'nS' 1 V-^) 

The most general solution to these n constraints (|5.27|) is given by 


n. 


(5.27) 


2 ( 71 - 1 ) 


9n 


II ^ = W/ '( X ) 2 - /(®)> 


(5.28) 


»=o 


where f(x) is the most general even polynomial of order 2 n — 2 


71—1 


/( x ) = ^bix 21 . 


(5.29) 


z=o 


Note that we have now recovered the solution to the classical loop equation that 
we obtained in section [|. We now repeat the procedure outlined there and define the 
Riemann surface E given by 


y 2 = W\x) 2 -f{x). 


(5.30) 
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The filling fractions S) then become period integrals of the meromorphic 1-form ydx 
over the 1 -cycle Ai of £ that encircles the i th branch cut 


Si = 


ydx 
. 2 ni 


(5.31) 


We can then argue that the change in the free energy due to bringing an eigenvalue 
from infinity to the i th cut is 


df 0 

dSi 


= / ydx. 


(5.32) 


'Bi 


Note here that the B cycles are non compact, so for (|5.32|) to make sense we have to 
introduce an ultraviolet cut-off A 0 in the integral which has been identified with the 
bare coupling of the gauge theory [Q. We comment here that there are only serni- 
classical arguments for equation (|5.32|) , and we have been unable to rigorously prove 
this as a consequence of the loop equations and 


5.3.2 


Contributions 


Once we have the form of the solution for uiq{x), we can substitute it in the loop 
equation, which is now a linear inhomogenous integral equation for uq(x), 


c 


dx' W' (x')uji{x') 1 

—- - -= u q {x)u-l[x) - — UJ 0 {x). 

2 m x — x' 2x 


(5.33) 


We can get a natural ansatz for uj\ from the string theory expectation that T\ 
should be a derivative with respect to So of jr 0) 

dBo 


Ti = q 


dSo ’ 


(5.34) 


where q is some constant which has to be determined. Inserting this into ( 5.12 ), we 
get 


d 2 j d dT 3 

gs~ o 


dV 
d 


S 


=q w 0 C o(x) ~ * 

duJo 2 q 

=q- 


(5.35) 


dSo x 

It is easy to see that solves the homogeneous part of the loop equation. The 
inhomogenous part of the loop equation is solved by — — if q = — 
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More generally, in the case of multi cut solutions, we could have added any solution 
to the homogeneous loop equations. This amounts to taking 

x = < 5 ' 36 ) 

i 

such that Qi = ~ |> However, corrections of the form for i > 0 should not 
be generated since these cuts represent U(Ni ) gauge physics. We will give a short 
perturbative discussion of this in the next section. 

We can extend this result to a single cut model with an arbitrary polynomial po¬ 
tential. We will use this to solve the SO/Sp Af = 1* theories. 


5.4 Counting Feynman diagrams with S 2 and MP 2 topology 


For a perturbative check of the relation 


T\ = ±q 


dfo 
dS n 


(5.37) 


we need to enumerate “ribbon” graphs in the’t Hooft (genus) expansion of the matrix 
model. Recall that the genus expansion is ordered by diagram topology, with diagrams 
of genus g and c cross-caps contributing at order g~ x = g^ 2+29+c - The coefficient q is 
related to the relative contribution of the planar (genus 0) diagrams which dominate 
at large M with the leading jj correction coming from diagrams with topology MP 2 . 

It is known that SO(2M) and Sp(M) matrix models are related by analytic contin- 

£J). Therefore, at 


uation M » — M (for the analogous gauge theory results see |]69| , 
even orders in the genus expansion, the contribution to the matrix model free energy is 
the same for both theories, while at odd orders the Sp(Ad) diagrams contribute to the 
free energy with an additional minus sign relative to SO(2M). This fact determines 
the sign in Q5.37p . 

Recall that 

X = v — p + l (5.38) 

where v is the number of vertices in the ribbon graph, p is the number of propa¬ 
gators and l the number of boundary loops. The Feynman rules are summarized in 
appendix |Q. Let us evaluate the first-order quartic diagrams in fig. [I]. The planar 
diagram has the value 


2 x (—Ym 3 

1! Ag s V2 mJ 
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(5.39) 







whereas the 


diagram with one twisted propagator contributes 

2 


-4 X iiL 

1 ! Ag s \2mJ 


M 


and the the MP 2 diagram with two twisted propagators contributes 

2 


1 X iiL (TXV 

1 ! 4 g s \2 mJ 

Using the fact, that S = %M, this shows that 


M 


Ti = ~ 


19Pn 


(5.40) 


(5.41) 


(5.42) 


4 

at the hrst order. 

We have calculated the Feynman diagrams for several higher orders and higher 
vertices and confirmed this relationship in those casesQ It would be nice to have a 
purely combinatorial proof of (|5.42|) , that would not rely on the loop equation. 

In order to describe a multi-cut matrix model, we need to use ghosts [[d 7|] to expand 
around the classical vacuum. In this prescription, one can think of the matrix model 
as several matrix models, which are coupled by bifundamental ghosts. Only one of 
those matrix models is actually an SO(Mq) / Sp(A4o/2) matrix model, the other matrix 
models are U(Mi) matrix models. The ghosts do not have twisted propagators either, 
so the leading contribution from the SO(M 0 ) / Sp(M 0 / 2) matrix model is again the 
same as for a single cut model. The loop equations still hold for the multi-cut model 
and we can extend the result to all orders. 


6 Computation of Effective Superpotentials 

In this section we combine the results of the previous sections to compute the effective 
superpotential of the dual gauge theories. We will find that it is necessary to refine the 
formula for the unoriented string contribution to the effective superpotential of [Q. 


6.1 Non-Perturbative Sector 


As discussed in [71,3], there is a non-perturbative contribution to the free energy which 
arises from the Gaussian integral: 


T nv = -dun G log ^ - log vol(G). 
_ 2 m 


( 6 . 1 ) 


“This relationship was apparently not known to mathematicians. 
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Planar diagrams: 








1 

Figure 1: S 2 and MP 2 diagrams with one quartic vertex, written in terms 
of twisted and untwisted propagators and as diagrams on MP 2 to show their 
planarity. Propagators that pass through the cross-cap become twisted. 
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In appendix B, following [|7TJ, we have included the large M expansion of the logarithm 
of the volume of the SO/Sp groups. We find that, for SO(M) when M is even, 


jrnp = l_jrnp + ]_jrnp + 
9 s 9s 


9s L 


s 2 log 


2ttS 


m 


S' 2 - + log 7T 


( 6 . 2 ) 


1 

9s 


S', 2 nS , 

-log- + - 1 + logTT 

2 m 2 


log 4) 


+ 


with a similar expression for M odd or G = Sp(M). We see that 


TT = T 


1 d^ p 
4 ds 


± 2 lo g 2 > 


(6.3) 


where the first — / + sign is for SO/Sp respectively. This is almost the same relationship 
as we found for the perturbative contributions (|5.42|) , but it is spoiled by the log 2 term. 
We have traced this term through the volume computation outlined in |7I| and found 
that it could be removed by a different choice for the measure on the maximal torus of 
the Lie group. 

It is the non-perturbative sector, specifically the coefficient of the S 2 log S term, 
that determines the number of gauge theory vacua, which is a main consistency test of 
the translation between matrix model quantities and the effective superpotential of the 
gauge theory. The number of vacua of a supersymmetric gauge theory is equal to the 
dual Coxeter number h of the gauge group m |73j . Therefore the total superpotential 
should lead to the conclusion that S h is single-valued. 

Open string physics tells us that the sphere contribution to the effective superpo¬ 
tential should be proportional to Qd 5 , the total charge of D5-branes, while the MP 2 
contribution should be proportional to Q 05 , the total charge of 05-planes. We can 
express this by refining the suggestion of [HI: 


37- 

W eS = Qd 5 + Q 05 Go ~ 27ri r S, 
We assume that Go is proportional to the total MP 2 free energy, 

Go = a(F7 + F p ). 


(6.4) 


(6.5) 


Proceeding with this result, we find that 


Weff = 


N a 

Y ± 4 


SlogS — -r S + ■ ■ 


( 6 . 6 ) 
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where the +/— is for SO/Sp respectively. Consistency with both the closed string 
result ( |3.10| ) and the gauge theory(] requires that we must have a = =[=4. Very re¬ 
cently 1 35j produced this factor |a| = 4 from a perturbative argument along the lines 
of |28|]. It was found to be related to the measure on the moduli space of Schwinger 
parameters, a quantity that is intrinsic to the gauge theory. Presumably, there is a 
similar explanation of this correction within the holomorphic Chern-Simons theory. 


6.2 The A/" = 1* Theories 

Following we can also consider the AT = 1* theories. Table [l] contains the results 
that are needed to write the partition function in the eigenvalue basis. In the large M 
limit, the discussion will entirely parallel that of j3j. Inclusion of MP 2 contributions to 
the superpotential for SO(2N) gauge group yields 

2N -2 1 

Weff=—g— nfi ( r ) - 2 T o u a(t). (6.7) 

This has extrema at 

= T o + k jfc = o ... 2JV - 1, (6.8) 

27V —2’ ’ v ; 

at which points the superpotential takes the critical values 

hF ( fe) ~ E 2 ((to + k)/{2N — 2)), (6.9) 


m 


complete agreement with ( |2.16|) . 


7 Discussion 

In this paper we have outlined a general scheme for computing the subleading con¬ 
tributions to the gauge theory effective superpotential from MP 2 diagrams in the dual 
matrix model of Dijkgraaf and Vafa. The methods involve an application of the higher- 
genus loop equations to determine the MP 2 correction to the resolvent, which allows 
us to compute the MP 2 contribution to the free energy of the matrix model. We then 
established a refinement of Dijkgraaf and Vafa’s relationship between matrix model 

7 Note that, after including a = y4, the effective superpotential naively suggests that for gauge 
group Sp(N/ 2), S N+2 is single-valued, whereas h = N/2 + 1. The resolution to this puzzle was 
explained in |6(i| l. Namely the Dl-string wrapped on P 1 has instanton number two in Sp(N/2). 
Properly accounting for this reproduces the Z 2 h chiral symmetry of the dual gauge theory. 
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quantities and the gauge theory effective superpotential which was necessary to obtain 
consistent field theoretic results. The computation of provides a gauge theoretic 
explanation of our prescription. 

There are many future directions that can be pursued. First, while we were inter¬ 
ested in an order g s (equivalently 1 /M) correction from unoriented diagrams, in the 
matrix model duals to gauge theory with matter in general representations, such as 
quarks in the fundamental, there will be 1/M correction^ arising from worldsheets with 
a single boundary. Our application of the loop equations should apply to this case and 


it would be interesting to use this to make contact with the results of 113,14,15,16,12 


Similarly, it would be of interest to examine higher-order corrections in the genus 
expansion, which have an interpretation as gravitational corrections si, mm 


it 


would be interesting to make contact between the loop equations and the Kodaira- 
Spencer equations of f74j], which also relate higher genus results to those at lower 
genus. It seems reasonable that there are higher genus forms of our relation between 
oriented and unoriented contributions at a given genus like 


F g , l oc 


dT a , o 
dS 0 ' 


(7.1) 


It would also be useful to obtain a deeper understanding of relations like ( |7.1| ) from 
the diagrammatic combinatorics. For example a diagrammatic proof of our conjecture 
seems to be possible. 
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A Matrix Integral Measures and Determinants 


In this section we collect some results on the group measure and adjoint action which 
are needed to do computations in the matrix models. 

A.l The Group Measure for General Matrices 

We wish to compute the Jacobian for the transformation from certain matrices $ to 
their eigenvalues. This can be derived by a group-theoretic argument. In terms of 
the Cartan generators Hi and ladder operators E a: for the algebra of the group G, 
satisfying 

[H i ,E a \=a i E a , (A.l) 

we can diagonalize a matrix $ 

$ = WAU, 

A = J>iT. < A - 2 > 


We will define parameters t a so that 


dU = 


^ ^ dt a E c 


U t* = —t 

o, L a L-c 


The infinitesimal variation of can then be written as 


=f/ t 
=t/t 

We now calculate the metric on the Lie algebra 


U 


dA + J2 dt a [A, E° 

a 

dA + dt a ( Acd ) E° 


U. 


Tr d& = £"? + £ dt a dtp £ x - a ‘ £ ^ E ° E t 

i a,(3 \ i /Vi / 


Using the identity 


(A.3) 


(A.4) 


(A.5) 


Tiy E a E? = C(r)S Q+ p fi (A.6) 

where C'(r) is a representation dependent constant, we can simplify the second term 
in equation (|A.5|) to 


C(r)£ 


£«'a, 


| dt 0 


(A.7) 











Up to numerical factors, the Jacobian is 


A(A) = ]; 


Q>0 




(A.8) 


A.2 The Induced Measure of the Af = l* Theory 


In this section, we calculate det (adj $ + im ) whose modulus squared appears in the 
calculation of the induced measure of the Af = 1* theory. In order to calculate the 
determinant, we go to a diagonal basis in which $ is an element of the Cartan subal¬ 
gebra ( |A.2|) . Then we solve the eigenvalue equation 


[<£>, A\+ iA = aA, 

where A is a completely general matrix in the Lie algebra 

A = ^h l H i + J2 t C *E a . 

i a 

We can compute 


[$, A]+iA — i ^2 hi Hi + S ** + i ) t a E a , 

i a \ i / 


so the eigenvectors and eigenvalues are 

A = H l , a = i, 

A = E a , a = (y^ j a l X i + i). 

i 

Up to numerical factors, the determinant is then 


det (adj $ + im) 


n E a* A i + im 1 

a>0 \ i ) 


(A.9) 


(A.10) 


(A.11) 


(A.12) 


(A.13) 


We list the expressions for the roots and the corresponding determinants for the 
different classical groups in Table [l]. 
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G 

Roots 

A A) 

det(adj A + i ) 

An -i 

n(Ai - -^j) 2 

i<j 

e i - e j (i 7 l j) 

n (A — A + ®) (A — A — 

i<j 

Bn 

Kj * 

±et ± ej (i ^ j), ±d 

n((Ai - a,) 2 +i)((a+a,) 2 +1) n(Af +1) 

i<j i 

c N 

n(A?-A?) 2 n\ ? 

Kj l 

73 (±ej ± ej) (i ^ j), ±V^e-i 

n((Ai - A,) 2 + 2)((A,; + \j ) 2 + 2) n(2A? + 1) 

Kj i 

D n 

nw-w 

Kj 

±e t ± ej (i 7^ j) 

EI ((Ai Aj) 2 + l)((Aj + Aj) 2 + 1) 

Kj 


Tabic 1: The roots and the formulae for the Jacobians and determinants of the adjoint 
actions for the classical groups. 

B Asymptotic expansion of the gauge group volumes 


We now compute the asymptotic expansion of the volume of the gauge groups, which 
normalizes the partition function of the matrix model and provides the nonperturbative 
contribution to the free energy. The volumes are given by [71]: 

2 n+1 (2tt) n2+n ~1 


vol(SO(2N + 1)) = 
vol(SO(27V)) = 


(27V- 1)1(2#-3)! 
V2(2tt) n2 


. 3!1!’ 


vol(Sp(27V)) = 


(27V — 3)!(27V — 5)!... 3!l!(7V — 1)! ’ 

2-w( 27 r)iv 2 + iv 


(B.l) 


(27V — 1)!(27V — 3)!... 3! 1! 

We are interested in the large TV asymptotic expansion of the logarithm of the 
volumes in order to compute the non-perturbative contribution to the free energy. 


Following [71], we introduce the Barnes function 


G 2 (z + 1) = T (z)G 2 (z), G 2 ( 1) = 1. (B.2) 

Using the doubling formula for T(^), 

r(2z) = 2 2z_1 7r _ ^r(^)r(^ + i), (b.3) 
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and ( |B.2|) , can evaluate the denominator of the volume factors 

G d (N) = (21V - 1)!... 3!1! = ^-^2 n ^G 2 (N + 1 )G 2 (N + |) (B.4) 

Using the Binet integral formula 

1 1 r°° 

log r(z) = (z- -)logz - z + - log 2vr + 2 J (B.5) 

the asymptotic expansion of G^C^) is 

N 2 1 3 1 

log 0 2 (1V + 1) = — log N - — log N--N 2 + -N log 2 tt + 0(1). (B.6) 

By expanding log(fV — a) for large N, we obtain 

log G d (N) =N 2 log IV + IV 2 (-^ + log 2) 

1 IN 

+ -NlogN - — log IV + — (\og4n- 1) + 0(1). 

Putting all of this together, we find that 

logvol(SO(21V + 1)) 

= — N 2 log IV + 1V 2 (^ + log7r) 

1 IN 

~ -NlogN + — logIV + —(1 + log 4 + logvr) + 0(1), 
logvol(SO(21V)) 

=-!V 2 logfV + fV 2 (^ + log7r) (B.8) 

1 IN 

+ -NlogN + — log IV + y(-l + log 4 - logvr) + 0(1), 

logvol(Sp(21V)) 

= —N 2 log N + 1V 2 (^ + log 7r) 

1 IN 

- -NlogN + — logiV + —(1 - lo g 4 + log7r) + 0(1). 

C Matrix model Feynman rules and enumeration of diagrams 

We want to perturbatively evaluate the matrix integral 

j d.* e£* , (C.l) 
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where the potential W is given by 


w($) = —■ 

3 = 1 ^ 

and $ is a real antisymmetric M x M matrix. We can write this as 


(C.2) 


exp 


-Tr 


m 

~2 


<l > 2 + \ 

^ 2 j 

3 =2 J 


(C.3) 


where m — g\. Expanding the exponential leads to traces of integrals of the form 


d<f>e 


—Tr 


2" w Cf> ... $ 


a 


a 


<9J„ 




d<f> exp 


' mini ^ u m^n^ 

This integral can now be evaluated, leading to 


—Tr —<f> 2 - — Tr J<f> 


j =o 


(C.4) 


2 ng s 


m 


M(M-l) 

2 


0 


3 


31 31 


— _9sL Tr /2 

e 8m 


j=o 


(C.5) 


Differentiating step by step gives rise to expressions like 


3 


( 9s j g S T J 2 


3J mn V2rn" mini '2m Jmknke 8m 


_ /r r C C \ 9 s J Qs j — jp^-Tr J 2 

^mn\ ^nm\ ) ^ ''m2^2 ' ' ' 277^ " m k n k ^ m 


JTs 

2m 

3s 


■ -/ ,, 


2m 


■ d 


2 ;/; 


i^mm.^nn^ Zirn^. ) ^ 




-J — J 

0 ^ rrm 0 ^ n 

zm zm 


9s 

2m 


d n 


-f^Tr J 2 

e Sm 


(C.6) 


The indices m, and rp are contracted in traces as given in the interaction which can 
be interpreted as forming vertices. The combinatorics can the be interpreted diagram- 
matically, that one must connect all the legs of the vertices in all possible ways with 
untwisted and twisted propagators. Each twisted propagator contributes a factor of 
(-!)• 

The rules for evaluating a diagram are then: 

• Each kind of vertex with multiplicity Vj contributes a factor of vtt (jj^~) Vj ■ 
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Each propagator contributes a factor of j 5 -. 


Each twisted propagator contributes a factor of (—1). 


Each index loop contributes a factor of M = 


25 

9s ' 


The combinatorial factor of a diagram can be computed by counting all topologically 
equivalent ways in which the legs of the vertices can be connected. This has some 
subtleties, since some diagrams with twisted propagators can actually be planar. To 


handle this, we make use of the technique described in [33] to draw unoriented diagrams 
(see also [[T5|,[76|] for recent work on non-orientable ribbon diagrams in the mathematical 
literature). 

An MP 2 can be drawn in the plane as a disc, where antipodal points on the boundary 
are identified. MP 2 diagrams can then be drawn on that disc with some propagators 
going through the cross-cap at the boundary. The propagators going through the 
cross-cap are twisted propagators, whereas all the others are untwisted propagators. 

We can now also draw a planar diagram on the MP 2 . If it has more than one 
vertex, we can push one or several vertices through the cross-cap without destroying 
the planarity, but all the propagators going through the cross-cap are now twisted 
propagators. This operation contributes a multiplicative factor of 2 V ~ 1 to the number 
of planar diagrams at each order v. See Figure [l] for the enumeration of diagrams with 
1 quartic vertex. 

Using the relation between p and the number of vertices v % of valency i according 
to 

1 ~ (C.7) 




IVi 


the contribution of planar diagrams to the free energy of the SU(M) matrix model is 
given by 


°o ,(n) oo i (n) 

T 0 = V = V U-(^)"(^)!™’M 2 - (1 -^, (C.8) 

J v\ K ng s y m ' v\ y ng s K m 


where the sum is over diagrams with v vertices of valence 2 n, d^ is the number of 
planar diagrams at each order, and l counts the number of boundary loops of the ribbon 
graph. The propagator for SU(M) theories is twice that of the SO/Sp theories. In 
the second line we have simplified using (5.38) and (C.7). The number of diagrams of 
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topology S 2 (i.e. planar diagrams) in SU(M) matrix theory with a quartic potential is 
given by 

rf(4) = (2t;-l)!12* 


= 2,36,1728,145152,.... 


(C.9) 


(v + 2)! 

We are not aware of explicit generating functions for other vertex valences 2 n, but 
these diagrams can be enumerated by computer to the desired order 0- 

If we now include twisted propagators (i.e. enumerate planar diagrams in the SO or 
Sp matrix models), there is an extra contribution to the set of planar diagrams coming 
from vertices that have been “flipped”, converting untwisted to twisted propagators 
according to the rule described above. 


T 0 = 



9n_\v (_9s _\\ nv , 
ng s v 2 m 7 


(C.10) 



1 (2v - 1)124" 

2 


2,72,6912,1161216,.... 


A similar expression exists for the MP 2 free energy 


?i = 




(C.ll) 


(C.12) 


Here the number of diagrams dfcf* is counted with a minus sign for each twisted prop- 
agatoif]. The relevant planar and MP 2 diagrams were enumerated by computer up to 
4 vertices with a quartic potential ITtree ~ 4* 4 , to 2 vertices with a sextic potential 
^tree 4> 6 , and for a single vertex with a potential of degree up to 16. The results 
are summarized in Table and verify the desired relation: 


iaPo 
2 dM ' 


(C.13) 


9 Gaussian Ensembles are matrix models that have been well-studied in the physics and mathematics 
literature. The Gaussian Orthogonal and Gaussian Symplectic Ensembles also contain non-oriented 
ribbon diagrams with twisted propagators, however the propagator is (T b °T^) ~ S ac Sbd + SadSbc > *-e-, 
there is no relative minus sign between the two terms. This corresponds to counting MP 2 diagrams 
with a positive sign always. Therefore the free energy of the Gaussian Ensembles differs from that of 
the Lie Algebra matrix models at sub-leading orders in the genus expansion. 
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Diagrams with quartic vertices: 


Gauge group 

Topology 

v — 1 

v = 2 

v = 3 

v = 4 

SU 

S 2 

2 M 3 

36 M 4 

1728M 5 

145152M 6 

SO/Sp 

S 2 

2 M 3 

72 M 4 

6912 M 5 

1161216M 6 

SO/Sp 

MP 2 

—3 M 2 

-144M 3 

-17280M 4 

—3483648 M 5 


Diagrams with sextic vertices: 


Gauge group 

Topology 

v = 1 

v = 2 

SU 

S 2 

5 M 4 

600M 5 

SO/Sp 

S 2 

5 M 4 

1200M 6 

SO/Sp 

MP 2 

—10 M 3 

-3600M 5 


Table 2: Contribution to the free energy of the SU/SO/Sp matrix models at 
planar and MP 2 level, for quartic and sextic potentials. The first few terms in the 
perturbative expansion are listed, corresponding to the number of diagrams with 
increasing number of vertices (equivalently loops). 
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